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ABSTRACT 

Precision Doppler measurements from the Keck/HIRES spectrometer reveal 
periodic Keplerian velocity variations in the stars HD 16141 and HD 46375. HD 
16141 (G5 IV) has a period of 75.8 d and a velocity amplitude of 11 m s -1 , 
yielding a companion having Msini = 0.22 Mjup and a semimajor axis, a = 
0.35 AU. HD 46375 (Kl IV/V) has a period of 3.024 d and a velocity amplitude 
of 35 m s -1 , yielding a companion with Msini =0.25 Mjup, a semimajor axis of a 
= 0.041 AU, and an eccentricity of 0.04 (consistent with zero). These companions 
contribute to the rising planet mass function toward lower masses. 

Subject headings: planetary systems - stars: individual (HD 16141, HD 46375) 

1. Introduction 

The two thousand nearest and brightest dwarf stars ranging in spectral type from late 
F through M are currently being surveyed by groups with a Doppler precision of ~10 m s _1 
or better. To date these surveys have resulted in the discoveries of 32 extrasolar planets 
(cf. Marcy, Cochran, and Mayor 2000, Marcy & Butler 2000, Vogt et al. 2000, Queloz et 
al. (2000), Udry et al. (2000), Noyes et al. 1997), including the first system of multiple 
planets (Butler et al. 1999) and the first detection of a transiting planet (Henry et al. 2000; 
Charbonneau et al. 2000). 

Remarkably, all 32 companions found from precision Doppler surveys have M sin i less 
than 8 Mjup, though companions with masses of 10-80 Mjup would have been much easier 
to detect. The companion mass function rises steeply toward smaller masses (Marcy & 
Butler 2000), but it turns over near 0.5 Mjup, presumably due to poor detectability. With 
conventional precision of 10 m s _1 , companions of 0.5 Mjup in 4-day orbits induce stellar 
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motion only a few times greater than such Doppler errors. The stars 51 Peg and HD 75289 
have the lowest known values of Msini , both 0.46 M JUP (Mayor & Queloz 1995, Udry 
et al. 2000). However with a precision of 3 m s -1 , one may explore the mass distribution 
below 1 M S at (= 0.298 M JUP ). If sub-Saturn mass companions occur less frequently than 
Jupiter-mass companions, the status of all Jupiter-mass companions as "planets" would be 
cast in doubt. Such a distribution of masses, peaked at ~1 Mjup, is not the case in our 
Solar System nor do theories of planet formation predict such a peak (eg. Lissauer 1995, 
Boss 1995, Levison et al. 1998). 

Here we report Doppler variations in HD 16141 and HD 46375 exhibiting Msini < 1 
M S at- Stellar characteristics and observations are discussed in the second section. Orbital 
solutions are presented in the third section, followed by discussion of results. 



2. Observations and Stellar Characteristics 

HD 16141 (79 Cet, HIP 12048, G5 IV) and HD 46375 (HIP 31246, Kl IV-V) are 
slowly rotating, chromospherically inactive stars, based on the weak Ca II H&K emission 
in their spectra. Their Hipparcos (Perryman et al. 1997) distances are 35.9 and 33.4 pc, 
respectively, giving them absolute magnitudes of My=4.05 and M v = 5.29, which places 
both stars 1.0 mag above the zero-age main sequence. The Hipparcos mission made 60 and 
63 photometric observations of HD 16141 and HD 46375, respectively, revealing that both 
stars are photometrically stable at the level of errors, ~0.01 mag. 

The mass of HD 16141 is M=1.01 M , and it has [Fe/H] = 0.02, as determined by 
high-resolution spectral synthesis (Fuhrmann 1998). The mass of HD 46375 is probably 
~1.0 M Q also. Its spectral type of Kl IV-V and color, B — V— 0.86, would normally 
indicate a mass of 0.9 M Q . But its metalicity appears to be high, [Fe/H]=~+0.34, as 
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judged from narrow band photometry (Apps, priv. comm. 2000). This high metalicity 
implies a higher mass than is associated with solar-metalicity stars of its color, giving it an 
estimated mass of 0.9 - 1.0 M . Here we adopt a mass of 1.0±0.1 M , for computation of 
Msini and semimajor axis. HD 46375 resides only a few arcminutes North- East from the 
bright nebulosity of the Rosette Nebula (which is over 1 kpc in the background), leaving 
some concern about contamination of the optical spectra and photometry of the star. 

The velocities of HD 16141 and HD 46375 have been monitored since 1996 and 1998 
respectively. The technique is the same as described in Vogt et al. (2000), using the Keck 1 
HIRES spectrometer (Vogt et al. 1994). The H&K lines near 3950 A provide a simultaneous 
chromospheric diagnostic (Saar et al. 1998), giving logi?'(HK) = —5.05 and —4.94 for HD 
16141 and HD 46375, respectively, typical for chromospherically quiet stars (Noyes et al. 
1984). The photospheric Doppler "jitter" of such stars is less than 3 m s _1 (Saar et al. 
1998). 

3. Orbital Solutions 

The 46 observations of HD 16141 are shown in Figure 1 and listed in Table 1. These 
observations have an RMS of 7.0 m s -1 , 2.5 times larger than the measurement error. A 
periodogram of these velocities (Figure 2a), reveals a dominant period at 75.6 days. The 
reality of this periodicity is supported by two tests. We broke the velocity set into its first 
half and second half, yielding separate periodograms with highest peaks at 75 d and 77 
d, respectively, both having a false alarm probability under 1%. Thus, both halves of the 
measurements reveal the 76-day periodicity. We also computed the false alarm probability 
of the 76 d periodicity in the full velocity set by using a Monte Carlo approach (Gilliland & 
Baliunas 1987). We generated 10 5 sets of artificial velocities drawn from a Gaussian error 
distribution while retaining the actual times of observation. None of these artificial velocity 
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sets yielded a periodogram peak as high as that actually found for HD 16141 (Fig 2a). 
Thus, the false alarm probability is less than lxlCT 5 . 

The best-fit Keplerian model to the velocities of HD 16141 is shown in Figure 1, and 
yields a period, P=75.82 d. The amplitude is, i^=10.8 m s -1 , and the eccentricity is e=0.28. 
The RMS of the residuals to the Keplerian fit is 3.2 m s _1 , similar to the median internal 
error of 2.8 m s" 1 . Adopting a stellar mass of 1 M Q , the companion has Msini =0.22 
Mjup, and the semimajor axis is 0.35 AU. A periodogram of the velocity residuals is shown 
in Figure 2b which shows no significant peaks, indicating that no additional companions 
are evident. 

The 24 observations of HD 46375 are listed in Table 2. Figure 3 shows observations 
obtained between 6 and 11 Feb 2000, revealing a three day period. A phased version of 
the entire data set spanning 515 d is shown in Figure 4. The same best-fit sinusoid is 
shown in both Figures 3 and 4. The RMS of the sinusoidal fit is 2.59 m s" 1 , slightly greater 
than the velocity uncertainty of 2.2 m s _1 . The best-fit Keplerian gives e=0.04±0.04 
(i.e. not significant), and yields an RMS of 2.44 m s -1 . The period is 3.024 d and the 
velocity semiamplitude K = 35 m s _1 . Assuming the host star is 1 M , the companion has 
Msini =0.25 Mjup and the semimajor axis is 0.041 AU. The orbital parameters for HD 
16141 and HD 46375 are listed in Table 3. 



4. Discussion 

HD 16141 exhibits the smallest velocity amplitude, K=10.8 m s -1 , reported for any 
planet candidate, thus warranting an assessment of stellar and instrumental sources of error 
Chromospherically inactive stars are intrinsically stable at the level of 3 m s _1 (Saar et al. 
1998). Further, the majority of our program stars exhibit no instrumental "drift" in the 
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velocity zero-point at the 2 m/s level during the last four years (Vogt et al. 2000). There 
is no plausible stellar time scale near 75 d except rotation. But spots are not important as 
HD 16141 is chromospherically quiet at Ca II H&K, and photometrically stable at 0.01 mag 
from Hipparcos. Thus, the most likely explanation for the Doppler variations is an orbiting 
planet. 

The eccentricity for the planet around HD 16141 is not yet well determined, 
e=0.28±0.15, and indeed cannot be reliably distinguished from circular without further 
measurements. However, all 21 planets orbiting beyond 0.2 AU (Marcy & Butler 2000) 
have eccentricities above 0.1, and thus this planet will constitute an interesting test case of 
eccentricities of a possibly low mass planet. 

The companion to HD 46375 is similar to the other "51 Peg-like" planets with their 
orbital periods of 3 to 5 days and circular orbits. Circular orbits are expected from tidal 
coupling with the primary star (Rasio & Ford 1996; Marcy et al. 1997; Ford et al. 1998). 
As the primary in the HD 46375 system has not been spun up, the mass of the planet is 
constrained to be less than ~15 M JUP (Ford et al. 1998, Marcy et al. 1997). The suspected 
high metalicity of HD 46375 of [Fe/H]=0.34 (Apps, priv. comm.) supports the suggestion 
that "51 Peg-like" planets are associated with high-metalicity stars (Gonzalez et al. 1999, 
Queloz et al. 2000) . 

The value of Msini =0.25 Mjup for the companion to HD 46375 is smaller than that 
of any known "51 Peg-like" extrasolar planet, with the smallest having been 51 Peg (0.46 
Mjup, Mayor & Queloz 1995) and HD 75289 (0.46 M JUP , Udry et al. 2000). This suggests 
that if orbital migration brings such planets inward (Lin et al. 1996), the process may 
continue to operate at masses near 1 M S at, pending knowledge of sini. 

With demonstrated precision of 3 m s _1 (Vogt et al. 2000), the Keck survey is currently 
capable of making 3 a detections of "51 Peg-like" planets down to the neptune-mass 
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range. Knowledge of the companion mass function of "51 Peg-like" planets down into this 
range will provide useful constraints on models that explain the formation and subsequent 
dynamics of "51 Peg-like" planets. 

The effective temperature of a planet in a 3.024 day orbit about HD 46375 would be 
~1400 K (Burrows et al. 1998). Henry (2000) reports that no transit occurs, implying that 
sin? < 0.992 . The expected amplitude of a transit signal is about 0.015 mags (Burrows et 
al. 1998, Henry et al. 2000, Charbonneau et al. 2000). 

The companions to HD 16141 and HD 46375 have the lowest values of Msini (0.22 
Mjup, 0.25 Mjup) found to date for extrasolar planets. The observed histogram of 
Msini shows a steep rise toward the lowest masses, consistent with a power-law, dN/dM 
oc M _1 (Marcy and Butler 2000). These new planets support the suggestion that the mass 
distribution continues rising to 1 M S at- Verification of any rise in the planetary mass 
function below 1 MsatwHI require more detections to account properly for incompleteness. 

We thank Kevin Apps for assessment of stellar characteristics. We acknowledge 
support by NASA grant NAG5-8299 and NSF grant AST95-20443 (to GWM), by NSF grant 
AST-9619418 and NASA grant NAG5-4445 (to SSV), travel support from the Carnegie 
Institution of Washington (to RPB), and by Sun Microsystems. This research has made use 
of the SIMBAD database, operated at CDS, Strasbourg, France. 
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Fig. 1. — Keck Doppler velocities, with error bars, for HD 16141. The solid line through 
the points shows the best Keplerian fit. The residuals having rms=3.24 m s _1 are shown at 
bottom (crosses) with arbitrary zero-point (horizontal line). 



Fig. 2. — a) Periodogram of the measured Doppler velocities of HD 16141. The dotted line 
shows the 1% false alarm level. The highest peak at 75.6 days has a false alarm probability 
less than lxl0~ 5 from Monte Carlo tests, b) Periodogram of the velocities of HD 16141, 
after removing the best fit Keplerian. No significant peaks remain. 



Fig. 3. — Keck Doppler velocities for HD 46375 from the 2000 Feb observing run. A 3 day 
periodicity is evident for this 6 night observing string. Measurement uncertainties are ~ 2.2 
m s -1 . 



Fig. 4. — Phased Doppler velocities for HD 46375. The same best-fit sinusoid is shown in 
Figures 4 and 5. The RMS to the sinusoidal fit is 2.59 m s -1 , essentially identical to the 
best-fit Keplerian. The period is 3.024 d and the semiamplitude is 35 m s -1 . Assuming the 
mass of HD 46375 is 1 M , the minimum (Msini ) mass of the companion is 0.25 Mjup, 
and the semimajor axis is 0.041 AU. 
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Table 1 — Continued 
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Table 1 — Continued 
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Table 2. Velocities for HD 46375 
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Table 2 — Continued 
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Table 3. Orbital Parameters 



Parameter 
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